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past  a synmetric  wedge.  Tulin's  (1964)  single-spiral  cavity 
model  is  used  to  represent  the  flow  and  conformal  mapping  and 
Riemann-Hilbert  technique  have  been  employed  to  obtain  the 
numerical  solution.  The  present  theory  shows  that  the  cavity 
length  for  the  gravity  solution  is  smaller  and  lies  entirely 
within  the  gravity-free  case.  Also,  the  change  in  the  drag 

coefficient  from  the  non-gravity  case  to  the  gravity  case  is  , 1 

small. 
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NOMENCLATURE 


A = the  onlv  non-zero  power  series  coefficient 
o 

Aj  = constant,  power  series  Jth  coefficient 

K2  = lergth  of  the  cavity  (fran  point  of  detachment  to  the 
rear  of  the  cavity  ) ( k2*=  C ) 

H(t)=  solution  of  the  homogeneous  Hilbert  problem 

**  •*-  - 

i.  = wedge  am  length 

PQ  ® reference  pressure  at  infinity 

p = pressure 

pc  = pressure  in  cavity 

q = magnitude  of  velocity 

qQ=  magnitude  of  velocity  at  infinity 

Q(t )=  solution  of  Riemann-Hilbert  problem 

t = auxiliary  half-plane 

t^=  real  t-plane  coordinate 

W = complex  potential 

x = horizontal  coordinate  in  the  physical  plane 
y = vertical  coordinate  in  the  physical  plane 
= semi  angle  of  wedge 
= normalized  conjugate  complex  velocity 
^ = dummy  variable 
^ = constant,  density  of  fluid 
<T  = cavitation  number 
X = dummy  variable 
cO  = log  ? 

<p  = velocity  potential 
^ = stream  function 
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1.  INTRODUCTION 


Only  a Tow  theories  Ivave  considered  the  Influence  of  a gravity 
field.  The  actual  magnitude  of  the  effect  of  the  gravity  depends 
upon  the  Proude  number  and  its  orientation  relative  to  the  flow 
direction.  In  general,  the  Froude  number  attained  In  water  tunnel 
tests  indicate  a small  gravitational  effects.  However,  the  effect  of 
the  gravity  may  become  Important  near  free  surface  and  certain  other 
applications.  The  purpose  of  this  work  is  to  study  the  transverse 
gravity  effect  on  a flow  past  a symmetric  wedge,  with  a finite 
trailing  cavity. 

Parkin  £ I }and  Street  [ have  presented  linearized  theories 
for  flows  in  a transverse  gravity  field  using  an  additional  approx- 
imation in  the  boundary  conditions.  Street  [»]  has  this  approximation 
to  be  quantitatively  equivalent  to  a first-order  solution  in  the 

p 

small  parameter  1/P  . In  all  these  theories  the  assumption  has  been 
made  that  the  ambient  pressure  is  constant  on  each  of  the  free  stream- 
lines, being  lower  than  the  reference  pressure  at  the  forebody  on  the 
upper,  and  higher  than  this  reference  pressure  on  the  lower  stream- 
line. 

The  case  of  a longitudinal  gravity  field  was  studied  by  Acosta^  4] 
and  Lenau  & Street  { 5 ] . Acosta  £ 4 } has  developed  a linearized  theory 
by  considering  the  flow  past  a synmetric  wedge  whereas  Lenau  &Street[s 
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have  developed  a non-linear  theory  by  representing  the  cavitating  flow 
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by  a modified  Riabouchinsky  model.  They  have  found  that  for  a gravity 
field  pointing  the  same  direction  as  the  flow  the  cavltv  is  shortened, 
while  for  a gravity  field  in  a direction  pointing  in  a opposite  to  that 
of  the  flow,  the  cavity  is  lengthened.  It  is  also  interesting  to  note 
in  these  theories  that  after  some  critical  values  of  the  field  strength, 
the  trailing  edge  of  the  cavity  becomes  cusped  rather  than  blunt.  As 
mentioned  by  Tulin  £ 6 j real  cavity  flow  involving  longitudinal  gravity 
fields  are  in  most  cases  unsteady  and  the  effect  of  the  unsteadiness  can 
easily  predominate. 

The  effect  of  gravity  on  the  stability  of  the  cavity  is  discussed 

by  Chang- Yi  Wang  f 7 J . 


2.  GENERAL  THEORY 


We  consider  the  steady,  tvro-dimensioned  irrotational  flow  of  an 
inviscid,  incompressible  fluid  past  a fully  cavitating  syimetric  wedge 
of  angle  2*  under  the  influence  of  transverse  gravity.  By  fully 
cavitating  flow  we  mean  a cavity  extending  along  the  entire  suction 
side  and  even  beyond  the  trailing  edge.  Bernoulli's  theorem  can  be 
written 

p + 1/2  fq''  + gfy  = constant  (1) 

where  p - is  the  pressure 

q - is  the  speed  of  the  fluid 

y - is  the  vertical  distance  between  a point  In  the 


! 


fluid  and  some  reference  elevation 
f -is  the  fluid  density 
g - is  the  acceleration  due  to  gravity 
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We  define  the  cavitation  number 


(?) 

(3) 


where  £ is  some  characteristic  length,  and  will  be  taken  as  the 
wedge  length  in  what  follows.  In  the  presence  of  gravity,  the  fluid 
speed  q can  be  expressed  in  terms  of  aQ  and  the  Froude  number  by 


J L. 

^ 1 i+r  i ' 


W 


We  also  define  the  following  symbols  : 
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a 

•d 


T 


z * x + i 3 

W = § + ‘ "p 


iW 


1 


u-iV 


(5) 
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The  method  of  solution  consists  of  two  basic  Darts  : 

(a)  Construction  of  conformal  mapping  between  the  physical  plane  and 
the  carrolex  potential  plane  and  a parametric  half  plane. 

(b)  Solution  in  this  plane  of  a well  posed  Riemann-Hilbert  problem. 

Towards  this  goal  we  can  determine  W as  a function  of  the  half- 
plane variable  t by  conformal  mapping.  Then  by  using  ttie  Riemann-Hi lbert 
technique  we  construct  urt)  explicitly  as  a function  of  t and  once  ) 
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is  known,  the  drag  coefficient  can  be  determined  using  Bernoulli's  theorem. 

The  general  solution  of  the  Riemann-Hl lbert  nroblem  in  the  upper- 
half  plane  is  well  known  [8  , ?].  Tn  this  ni-oblem  if  the  imaginary  part 
of  an  analytic  function  Q(t)  is  known  at  all  points  on  the  boundary, i.e. , 
on  the  real  t-axis,  then  Q(t)  is  obtained  in  the  entire  upper-half  plane 
from  the  formula  : 


J-oo  J=0 


-(7) 


where  Aj  are  real  constants  and  J represents  all  positive  integers. 

Now,  it  remains  to  relate  to  Q(t)  so  that  Im  Q(t)  is  known  at 
every  point  on  the  real  line. 

We  map  W-plane  onto  the  t-plane  so  that  the  wedge  is  mapped  onto  a 
finite  segjnent  of  the  real  axis,  and  the  cavity  boundary  onto  the 
remainder  of  the  real  axis.  Then,  we  have  boundary  conditions  along  the 
real  axis  in  terms  of  <oit)  . 

We  can  convert  these  conditions  to  those  of  the  Riemann-Hilbert 
problem  follwoing  the  general  method  of  solution  reported  by  Cheng  and 
Rott  0®1  The  original  function  ooCt)  is  alternately  purely  real  or  purely 
imaginary  on  the  real  axis,  hence  we  construct  an  auxiliary  function  H(t) 
which  makes  the  quotient  Q(t)=  known  on  the  entire  real  axis,  and 

Hi*/ 

always  imaginary  as  reauired.  It  is  to  be  noted  that  H(t)  is  a solution 
of  the  homogeneous  mixed  boundary  value  problem;  i.e.,  the  prescribed 
values  of  tJrtJ  are  all  zero  along  the  rea]  axis.  The  choice  of  the 
homogeneous  solution  H(t)  is  quite  arbitrary. 
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A choice  of  H(t)  satisfying  our  renuirement  is  ( see  Fig.  3 
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H(t)  = i [(t-t^t-t^J 


1/2 


(8) 


where  t-^  and  t correspond  to  the  ends  of  the  wedge. 
Let 


3m|(H(t)j]  >0 


on  the  real  axis  for  t<t^,  and  select  a branch  cut  on  the  real  axis 
interval  (t^,  t 2)  so  that  H(t)  is  single  valued.  Hence  £J«>  can  explicitly 
be  constructed  by  the  use  of  the  Riemann-Hi lhert  solution  where 
w£*)=  H(t)  Q(t)  . 


3.  ANALYSIS 

The  solution  to  our  problem  is  based  on  an  iterative  procedure 
suggested  by  Larock  and  Street  £ 1 1 } which  employs  the  non-gravity  solution 
as  the  initial  trial  solution  ( The  non-gravity  solution  Is  obtained  frcrn 
Ref.  12  ).  We  are  using  the  same  cavity  model  ( single  spiral  vortex 
model  ) as  Is  used  in  the  non-gravity  case.  There  are  two  distinct 
numerical  advantages  in  using  this  particular  model  as  discussed  in 
Ref .11  ( see  also  figures  1,2  and  3 ). 

Let  us  denote  yL  and  yn  be  the  respective  values  of  y(t)  for  the 
lower  and  upper  streamlines  beyond  whclh  y(t)  becomes  constant  (Fig  1 ) 
and  corresponding  values  of  t are  t^  and  ty  . Then,  for  t£t^  , 
y(t)=  yL  ; for  t >,  ty  , y(t)=  yy  . 

In  this  problem  the  gravity  acts  parallel  tc  the  y-axis  and  in  the 
negative  direction.  The  formulation  can  now  be  completed  by  starting 
the  boundary  conditions,  ""hey  are  : 


1.  q = qc(  1 - 

4* 


J *1_  )1/2 


»+«■  Fx  %■ 

l0g  * ~ * ^ C l+<r) 


on  the  cavity  and  as  z — » ®o 
Tm  ( ) = 0 , t=  ik 
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2.  The  combination  of  the  body  and  cavity  must  bo  closed.  The 
equivalent  statement  is  that  the  net  source  strength  must  be  zero  (j  2^ 

3.  Lastly,  the  flow  cannot  contain  non-integrable  singularities 
on  the  slit  or  have  multiple  values  off  the  slit  |j  3 J. 

Conditions  1-3  are  sufficient  to  determine  the  flow  field. 

Using  the  Ref .12  we  can  construct  the  following  boundary  value 
problem  : 

I 2 3 It) 


R.M-  + -fr2!52] 


-<>0  4 -L  <■  - 1 


dm  [ « I = - «*  . 


■ I 4 t < O 

o < -i  < I 
oo 


(9) 


where  F is  a prescribed  Froude  number;  , IT,  and  y(t)  are  respectively, 
the  semi -wedge  angle,  cavitation  number,  wedge  length  and  the  ordinate 
of  the  free  streamline  for  the  gravity  problem. 

In  order  to  apply  the  Riemann-Hilbert  technique,  we  now  form  [l2] 

as  follows  : 

-i/ 


Hit)  J 

= “z"  hO'iTr 

l 

F1 

! 2 

-t*>  < t 4 - | 

•Jr*  [ 

uUt)  T 

Hit)  J 

= l+t)J 

*4 

t 

- 1 < t 4 O 

[ 

U (t)  | 

Hi*)  J 

= -<*((  l-t)C/+t)J 

'4 

f 

0 4t4| 

(U  f 

aw  I 

= j-  8,4  r, . _l_  • 

21u) 

|r 

1 1 ( 1 ii \ <4  1 \ 

HH)  J 

2 0 1 i+r  f* 

i ■ 

| < t £ 00 

(10) 
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The  Riemann-Hilbert  power  series  will  still  give  only  one  term  as  in 


Ref .12  . Hence 


^ ~Oo  ) A ^ 

+ ( ° «K  f ^ 

11  7_,  IT  J0Ln.-i) 

+ _!  ("hit'-iTF  ?^’J  i.,  , , } 

j,  Cn-tXil1- .)‘/a  + J- -(11) 

Here,  two  terms  have  been  added  to  the  non-gravity  expression  for  vd) 
i.e.,  the  integrals  over  (- oo  , -1)  and  (1  , <»  ) . 

In  order  to  evaluate  the  extra  integrals  we  will  follow  the  Larock 
and  Street  method  [ll  J.  Recalling  that  y(t)=  yL  , for  t £ t^  ; and 
y(t)=  yy  , for  t t^j  . We  obtain 


M'-fSr  , . , , , 

L lvt,(.V-V  ^ fcjO—jr 


2?(n) 


+ f'W1-  drji**£J#( 


l co  cVtK’0-  O*  i i ) Hl 

L 


Likewise, 


ITT- 71  JL  ] . _ . I 1 zjfir  \ . 

(OO  j -a  r.  I * 7 

• ] + Aj  0 - i+c  p*  Jt  J M 

i\j  /|  Cl-t  )C^1- I ) 


The  integrals  over  (-00,  t, ) and  (t^  ,°°)  can  be  evaluated  exactly. 


The  unknowns  of  the  problem  are  the  same  as  In  the  non-gravity  case.  We 
can  rewrite  Eq.  ( 11  ) as 
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toi-t)  = H c-t  j f — L_  c i - — | u — *1 

[ 2TT  t0J  i*r  P 9.  ')  (n-t)U'.*)  *• 


^ in  >+<r  p*  * ,)* 

+ J*  . r»  M * f1  ^ 

Tf  J-|  C^-t)Ci-^*)y'i  It  Jo 


.L^a00  o(« 


'-( U>-^kV> 


iw  cwcv->*  2t; 


+ J_  + A 

cvtJCY-')1'* 


r"n  the  real  t-axis,  taking  the  proper  branch  of  H(t)  in  the  different 
region,  we  obtain  [l2  J 


°jLt)  r ~ ir['i7t  ' fy('!+U-  t‘)V*  )J  + AoO-t2)^  4 JAct) 

+ ^J^1'  1JF  pr  J r 6 U'+t/*0-tL)^_  ct-tK  i+-tL)^  | 

2 • -**>*  1 ,Jl(i*t)»i(,-tt)H+(l.0K0+tt)»4j 

~ 2*J  { t-  id-t*)14  jj 


+ " ^ 1 ?*  it) 

2 K (.  I - t 2 ) '/* 


J j_  r 


~r 

- £,<1  j ci+t)54 ( i-tu/z-  o-^ci+-tur  i 

<7  l Ci+t)^(|-  tu)^+ 


- I ^ -t  ^ o 
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a 

- + l It  + 2 i oYctom  c t2- 1 )/aJ  i L * + i-jct) 

t%0-,Tyfiru ) F i i c.^Vn^-n-oV+tQ4  i 
17C  t-t1-  I )**•  CH-t)*(l-tO*  +Cl-t)»(l+-tO>4  > 

- Uj  { t ± Ct'- . £ ) J 

+ 1'xi2  r p f crtt^Cl-to)'4*  7 

+ i.J  { -t  ± (t1-/)4}  ] 


for  -£  4 - I , -tr*k<  (+)  S^n 


where 


^ ^ *1  ■ 1 ^ -t+t  lH-t)(^.  |)l/* 


_ j 1 h<-'~T *? 

-tL  C^-i  )'^  -* 


F»r  -t>  i 


(17) 


^ ^-t)CV-0^  “ JtL  (ftHf-i)^ 


-u-t 


j_  ajffoh 
F»  X , 
c*t-  * > (t-l  ) v*- 


t < -i 

(18) 

As  in  the  non-gravity  case,  we  have  two  unknowns  Aq  and  K.  Since  boundary 
conditions  at  infinity  (t=ik)  still  hold,  the  following  integrals  (using 
the  notation  of  Larock  and  Street  ) will  be  found  useful  : 

Aj  ^ ’ H-r  p*  % j 


I|  (H1 . 1lx)  - | 


\ (Xl+  K1)  (V-0 


Ij.  , *u) 
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A numerical  Investigation  shows  that  the  integrals  (Eq.  12  ) over  the 

interval  (ty  ,oo  ) and  (-00  , t ) ( these  integrals  represent  the  spiralling 
region  of  free  streamlines  ) have  very  little  effect  on  the  other  integrals. 
Therefore,  in  further  calculation  we  will  neglected  these  integrals. 

The  two  equations  for  determining  the  unknowns  (similar  to  non- 
gravity case,  Ref.  1 2 ) are  : 


1+  K*  )Vx 


2 It 


{ I|  ('»  ±u)  - I,  (-fcL  ,-»)}  + A*(i+k*)' 


(19) 


re  K( K*>  1 ) Vx  + in.  1 T|  C',tu ^ -  H.  I,(-tL#-i)} 

f 2k(  j IZCI,  4a)-  I2  (-fcL;-i  ) ~ 


0 

(20) 


From  these  two  equations  we  can  always  find  A , K for  a given  oC  , tr,tT  and 

O JL# 


H.  CALCULATIONS  FOR  WEDGE  LENGTH,  DRAG  AND  LIFT 


The  expression  for  the  wedge  arm  length  is  : 


j \cU  \ = j 2KI(k\.) 
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Jiw 

p+O-t*)'* 
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The  expression  for  the  lift  and  drag  is  given  by  : 

•p  + iL  = - if  + i a 

A * 


l''rom  tho  Bem<>ull  1 'a  equat  ion  we  have 
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i+  (T  pi  JL  J 


refining  the  lift  and  drag  coefficients,  respectively,  as 


CL  = 


'4  X * 


CD  = 


We  have 


--- y*  U* 
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Let 


Then, 


TC  = 
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°L 


zK*(k%o  C JtCt)r  iW/%d* 

* -t,  ' J et‘+kl>z 


1 Sm  ^(l-l?F7I  S"d)  + ( lCj  - ICl)  **"  <*_ 

( IC|  + Ic»  ) * - - 


(21) 

(22) 


where  Ip  is  the  value  of  I_  with  limits  0 to  1;  and  I.  Is  the  value 

C ^ 2 

of  1^  with  limits  -1  to  0 . 

The  parametric  representation  of  the  free  streamline  Is  given  by 

* ~w£t)  dvj 

where  co(t)  is  given  by  Eq.  ( 14  ) and  tQ=  1 for  the  upper  streamline 


Z-  = j e " -42-  *fc 


(23) 


and  t = -1  for  the  lower  streamline, 
o 


5.  NUMERICAL  PROCEDURE 

The  Iterative  method  which  we  have  used  in  finding  the  numerical 
solution  to  the  present  problem  is  described  in  detail  by  Larock  and 
Street  £ 1 1 j . The  iterative  procedure  works  in  the  following  manner  : 

(1)  We  determine  all  the  parameters  and  the  physical  plane,  especially 
y(t)  on  the  streamlines,  from  the  non-gravity  case  ( Ref.  12  ) for  a 
given  K and  <S  . 

(2)  Using  the  input  data  from  step  (1)  where  appropriate  1^,  I2  and  fiQ  (Eqs.19, 

20)  are  then  computed.  Then  we  determine  J(t)  , y(t)  and  J^(t)  (Eqs.  18,23  & 15) 

(3)  Another  gravity  solution  is  obtained  using  the  data  from  the  step  (2) 
as  input  data. 
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(4)  Step  (3)  is  repeated  until  the  results  of  two  successive  integrations 
differ  by  less  tiian  sane  specified  amount;  at  this  point  we  stop  the 
calculations  and  gravity  solution  is  obtained. 

The  method  just  described  is  direct  and  simple  to  apply.  However, 
there  is  one  drawback  to  this  method,  namely,  that  it  is  very  time- 
consuming.  Even  computation  of  one  iteration  takes  a considerable  time 
on  a high  speed  computer. 

First,  y(t)  or  some  approximation  to  y(t)  is  assumed  to  be  known 
as  a continuous  function  of  t.  The  function  J(t)  and  Jg(t)  must  also 
be  determined  as  continuous  function  of  t since  they  occure  in  various 
integrands.  It  is  not  possible  to  compute  the  value  of  these  functions 
for  all  t because  it  would  be  a very  tedious  and  expensive  process . 
Therefore,  y(t),  J(t)  and  J|(t)  are  computed  for  a selected,  finite 
number  of  points  t^  . Then,  by  using  interpolation  methods  these  functions 
are  approximated  between  the  points  t^. 

In  numerical  calculations,  proper  care  must  be  taken  at  the  singul- 
arities of  the  integrands,  even  though  a singularity  is  integrable. 

In  Eq.  18  , valid  for  t < -1  , the  integrand  over  the  interval 

(1,  tjj)  has  a square  root  singularity  at  1.  We  treat  the  remainder  of  the 
integrand  as  a constant  over  the  interval  (1,  l+£  ) and  integrate  the 
singular  portion  : 
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For  a specified  t and  £ , the  term  can  now  be  computed  numerically. 

So  long  as  0<  Z < 1,  the  specific  choice  of  fc.  is  not  importnat. 

Likewise,  we  can  treat  the  integral  over  the  interval  (t  , -1) . 

L 

As  pointed  out  in  Ref.  12,  the  points  where  the  streamlines 
are  vertical  at  the  end  of  the  cavity  can  be  found  as  a function 
of  t by  finding  the  first  zero  of  the  integrand  of  equation  (4.1*0 
in  Ref.  12.  We  call  this  point  ty.  Since  we  are  considering  a 
symmetric  wedge  flow  in  the  non-gravity  case,  the  value  of  tL 
will  be  the  negative  of  ty.  ( The  roots  which  make  the  integrand 
zero  are  obtained  by  a bi-section  procedure.  ) 
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6.  DISCUSSION 

Figures  4 and  5 show  tho  effect  of  the  Froude  number  on  the 
cavity  shape.  The  cavity  for  the  gravity  solution  is  smaller  and 

lies  entirely  within  the  gravity-free  case.  Also,  the  change  In 
the  drag  coefficient  from  the  non-gravity  case  to  the  gravity  case 
is  small. 

From  the  present  theory  there  results  a negative  lift 
coefficient.  Intuitively,  it  seems  that  it  will  become  positive 
only  when  the  cavity  is  so  short  that  its  downward  pulling  effect 
is  balanced  by  the  buoyant  force  of  the  free  body. 

As  discussed  by  Larock  and  Street  [")  and  Parkin  £ I J it  can 
be  shown  in  the  present  case  also  that  the  effect  of  the  gravity 
Increases  with  decreasing  cavitation  number.  ( This  was  not  done 
because  of  computer  time  constraint  . ) This  simnly  means  that 
longer  cavities  are  more  affected  by  the  gravity  than  the  short  ones. 
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